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Abstract 

^ Like any fluid heated from below, the atmosphere is subject to ver- 

l/-) tical instability which triggers convection. Convection occurs on small 

time and space scales, which makes it a challenging feature to include 
in climate models. Usually sub-grid parameterizations are required. 
Here, we develop an alternative view based on a global thermody- 
namic variational principle. We compute convective flux profiles and 
temperature profiles at steady-state in an implicit way, by maximizing 
the associated entropy production rate. Two settings are examined, 
corresponding respectively to the idealized case of a gray atmosphere, 
and a realistic case based on a Net Exchange Formulation radiative 
scheme. In the second case, we are also able to discuss the effect 
of variations of the atmospheric composition, like a doubling of the 
carbon dioxide concentration. 
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1 Introduction 



Differential heating in the atmosphere drives atmospheric motion of primary 
importance for climate. On the horizontal direction, the meridional transport 
of energy effected by the atmosphere (and the ocean) moderates the tempera- 
ture gradient as compared to the conditions which would prevail if the planet 
were in local radiative equilibrium. On the vertical dimension, differential 
heating may trigger convection, which affects the lapse rate on the long term. 
Since convection occurs sporadically and on relatively small spatial scales, it 
is necessary to parameterize its effect in global climate models. As this is 
an important source of error, it would be desirable to investigate methods to 
compute the mean convective fluxes with the least possible adjustable coeffi- 
cients. In fact, horizontal heat transport suffers from the same problem: due 
to the turbulent nature of the transport, subgrid-models are also required. In 
this context, it was suggested by Paltridge (1975) that one may use a global 
thermodynamic variational principle to compute the atmospheric heat trans- 
port. More precisely, he conjectured that the meridional fluxes of energy may 
arrange so as to maximize a function, which corresponds to the rate of en- 
tropy produced by the meridional transport, subject to a global constraint. 
This hypothesis was criticized, rightly, for a number of reasons. First and 
foremost, there is no theoretical basis for the maximum entropy production 
(MEP) principle — see Ozawa et al. (2003); Kleidon and Lorenz (2005); Mar- 
tyushev and Seleznev (2006); Dewar et al. (2012) for general introductions. 
Some attempts were made to provide one, notably by Dewar (2003); Dewar 
(2005), but as of today a convincing proof is still lacking (see for instance the 
comments by Grinstein and Linsker (2007)), although connections exist with 
the phenomenological, well-established, theory of non-equilibrium thermody- 
namics (Groot and Mazur 1962[ ), and in particular with the Onsager theory 
which applies close to equilibrium (Onsager 1931). In addition, the Paltridge 
model suffers from a number of ad-hoc hypothesis. In spite of this, his results 
appear to be relatively robust (Grassl 1981; Gerard, Delcourt, and Francois 



1990 O'Brien and Stephens 1995 Herbert et al. 2011). In the case of the 



meridional heat transport, even though it may seem counter- intuitive to be 
able to reproduce a reasonable temperature distribution with no explicit rep- 
resentation of the dynamics of the atmosphere and ocean, one may argue that 
finding the correct shape for the meridional transport is not so surprising as it 
can be proved by other means that this feature is constrained by global ther- 
modynamic properties (Stone 1978). On the other hand, we know that heat 
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transport on the vertical involves different physical processes, characterized 
by faster time scales and smaller spatial scales. In particular, the hydrological 
cycle (Pauluis and Held |2002b| Pauluis and Held |2002a| Pauluis |2011[ ) and 
baroclinicity (Lucarini, Fraedrich, and Ragone 2011) play an important part 
in the thermodynamics of climate on the vertical dimension. Nevertheless, 
it seems appealing to try to apply the maximum entropy production princi- 
ple to compute vertical heat fluxes without making use of explicit sub-grid 
parameterizations. Along these lines, Ozawa and Ohmura (1997) obtained 
realistic temperature and convective flux profiles using a very simple model 
for the vertical structure of the atmosphere. In particular, they assumed 
that the atmosphere absorbed radiation as a gray body. The advantage of 
this unrealistic assumption is that radiative computations are straightfor- 
ward. These results were essentially confirmed by further indirect studies in 
similarly idealized frameworks of a gray (Pujol and Fort 2002) or semigray 
(Pujol 2003) atmosphere. Hitherto, the only study considering the full, non- 
gray character of atmospheric radiation focused on the role of water-vapor 
and clouds, treated with the MEP principle, while the thermal structure of 
the atmosphere was modeled with standard energy-balance methods (Wang, 
Nakajima, and Shi 2008). On the other hand, sensitivity experiments in a 
General Circulation Model (GCM) with respect to coefficients involved in 
sub-grid parameterizations, in which the entropy production rate was con- 
sidered as a diagnostic quantity, did not support these results (Pascale et al. 



2011). Similarly, the results of a 2D zonal-mean variational problem in which 



the diabatic heating field is extracted from a GCM show less agreement on 
the vertical dimension than on the horizontal (Pascale et al. 2012). In this 
study, we propose to compute MEP temperature and convective flux profiles 
for a realistic atmosphere, in the context of a predictive model formulated as 
a variational problem. To this end, we have developed a new radiative model, 
based on the Net Exchange Formulation for radiation (Dufresne et al. 2005). 
This model is a narrow band model based on the Goody statistical model, 
taking into account absorption by water vapor and carbon dioxide. We have 
solved numerically the variational problem for various standard water-vapor 
profiles. These solutions constitute the first MEP vertical profiles obtained 
with a realistic radiative scheme. In particular, they do not exhibit the 
vertical instability that appeared in previous studies with gray atmospheres 
(Ozawa and Ohmura 1997). The model also allows for further investigations, 
like the behaviour of the vertical profiles when the radiative properties of the 
atmosphere vary. Specific attention is devoted here to the case of a variation 
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in the carbon dioxide concentration. 

The outline of the paper is as follows: after briefly presenting MEP pro- 
files in the case of a gray atmosphere in section |2j we show in section [3] 
the profiles obtained with a new, realistic radiative scheme, based on the 
Net Exchange Formulation (NEF). In section |1J preliminary results about 
the climate sensitivity of the NEF-based MEP model are shown, while the 
conclusions are presented in section |5j 



2 MEP Profiles for a gray atmosphere 

Let us first briefly treat the idealized case of a gray atmosphere, following 



Ozawa and Ohmura (1997). We consider an atmospheric column made of 
TV layers. Each box has a temperature Tj. At the interface between layer 
i and i + 1, we note Fg(i) the net incoming (downward) solar flux, F L {{) 
the net outgoing (upward) longwave radiation, and Fc(i) the net (upward) 
convective flux. At steady-state we must have at each interface: 

F s (i) = F L (i) + F c (i). (1) 



2.1 A gray radiative model 

The radiative fluxes F$ and Fl depend on the composition of the atmosphere 
and on the temperature profile. 

2.1.1 Shortwave Radiation 

Let us assume here that the atmosphere is purely absorbing in the shortwave 
domain (no diffusion) and that the absorption is continuous, so that the solar 
fluxes satisfy the Beer-Lambert law: 

dF s (r) = -F s (r)dr, (2) 

where r denotes the (shortwave) optical depth. We choose atmospheric levels 
of equal shortwave optical depth, with total optical depth t$, so that 

F s (i) = F s (iV)e-( 1 -^) Ts . (3) 

The top-of-the-atmosphere net solar flux is given by Fs(N) = (1 — a p )S/A, 
where S = 1368 W.m -2 is the solar constant and a p = 0.3 the planetary 
albedo. The total shortwave optical depth is chosen so that the net surface 
solar flux is F s (0) = 142 W.m~ 2 , which gives Ts ~ 0.524812. 
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2.1.2 Longwave Radiation 



For the longwave fluxes, we also assume that the atmosphere is purely absorb- 
ing (no diffusion), with a linear extinction, in the framework of the Eddington 
approximation (two stream aproximation (Liou 2002)). The radiative trans- 
fer equation reads 



r/T 4 

4,— = 3F lW - 



d 2 F L (r) 
dr 2 



(4) 



where r denotes the longwave optical depth. For atmospheric levels with 
equal longwave optical depth, the radiative transfer equation can be dis- 
cretized as 



N 



4r L 



NJ 



Fr(i) - FrXi 



Ft a 



(5) 



where t l is the total longwave optical depth, a the Stefan-Boltzmann con- 
stant, and 1 < i < N — 1. Besides, 



oTt - aT 4 



<TT; 



N 



F L (N) N 
2 4r L 



(F L {N - 1) - F L {N)) . 



(6) 
(7) 



These equations can be recast in a simple matrix form. Let a = N/2tl. 
We define the iV + 1 dimensional square matrices A and S and the N + 1 
dimensional vectors T and F: 
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The radiative transfer equations now simply read 



AT = SF. (10) 

The matrices A and S being clearly invertible, one may at leisure express T 
as a function of F and conversely. 

We further assume that the shortwave and longwave optical depths are 
proportional, so that we can choose atmospheric levels which have at the 
same time equal shortwave optical depth and equal longwave optical depth. 
Then Eqs. [3] and 10 hold simultaneously. 



2.2 The MEP variational problem for convection 

The thermodynamic equation reads 
dT 

c p -^ = V.(F s -F L -F c ), (11) 

or in the discrete form: 
BT 

c v -^ = F s (i) - F s (z - 1) + F L (z - 1) - F L (i) + F c (z - 1) - F c {i). (12) 

At steady-state, this is equivalent to F s {i) = F L {i) + F c {i) for 1 < i < N. 
In these equations, the convective flux Fc is not determined. One possi- 
bility would be to build a sub-grid parameterization to compute it. Al- 
ternatively, we can formulate the Maximum Entropy Production variational 
problem, which states that, at steady-state, we should choose the values 
of Fc(l), ■ ■ ■ , Fc(N) maximizing the entropy production associated to con- 
vection. The (thermodynamic) entropy production associated to convection 
reads 



V F n (i - I) - Fn(% 



a = 

iv 



(14) 



i=0 



T 
1 



Ti+\ Ti 



where Fc{— 1) = 0. Mathematically, a can be seen either as a function of the 
A^+l variables F c (0), . . . , F C (N), because the temperatures can be expressed 
as T = A _1 SF and the longwave fluxes F are themselves functions of the 
convective fluxes at steady-state through the relation Fl{i) = Fg(i) — Fc(i), 
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Table 1: Gray atmosphere model parameters. 



S a p t s t l 

1368 W.m" 2 0.3 0.524812 2 -4 



or as a function of the N + 1 variables T ,Ti,. . . , T/v by reversing these 
relations. Let us for instance see a as a function of the convective fluxes, so 
that the MEP variational problem simply reads: 

max {a(F c (0),...,F c (N)),Fc(N) = 0}, (15) 

(F C (0),...,F C (N)) 

where the global constraint Fq{N) = enforces the steady-state condition. 
This is a constrained variational problem, but in practice it can be solved 
easily by deleting the term corresponding to % = N in the expression for the 



entropy production (14) and in the F vector. Alternatively, considering a as 



a function of the temperatures, the variational problem reads 

max a(T ,...,T N ). (16) 

(T ,...,T N ) 

This is an unconstrained variational principle (the constraint Fc(N) = can 
again be enforced simply by deleting the term corresponding to i = N in the 



expression for the entropy production (14)), and as such it is easy to solve 



numerically. Contrary to the two-dimensional horizontal case (Herbert et al. 



2011), the treatment of the steady-state constraint is thus straightforward 
here. The reason for this is that the formulation of the problem as a one- 
dimensional model with fluxes at the cell boundaries reduces the condition 
- vanishing integral of the divergence — to a condition on a single variable. 

2.3 Results for a gray atmosphere 

Given a set of numerical values for the parameters of the models, solving 



the model reduces to solving the variational problem (16). The parameters 
of the model are the solar constant S, the planetary albedo a p , the short- 
wave optical depth r$ and the longwave optical depth tl- Their values are 
recorded in Table [T] The variational problem can be solved using standard 
optimization routines. A unique maximum is found. The resulting profiles 
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Figure 1: Left: Temperature profiles at MEP state (solid lines) and at radia- 
tive equilibrium (dashed lines) for various longwave optical depths: Tl = 2 
(blue), tl = 3 (red) and tl = 4 (yellow). Right: Convective flux profiles at 
MEP state (solid lines) and for isothermal profiles (dashed lines) for various 
longwave optical depths: tl = 2 (blue), tl = 3 (red) and tl = 4 (yellow). 



are shown on Fig. [I] for various choices for the longwave optical depth. Note 
that the convective fluxes slightly differ from those obtained by Ozawa and 
Ohmura (1997) with the same model: we find an inversion of the convec- 
tive flux profile in the lower troposphere while they had found a constant 
profile in this region. This profile does not appear as a maximum entropy 
production state in our model, not even a local one. A possible interpreta- 
tion of this behavior is that the entropy production rate is locally negative 
in the inversion zone for the unconstrained variational problem. This locally 
negative entropy production rate is compensated for by a larger (positive) 
entropy production rate in the overlying regions, resulting in an overall en- 
tropy production rate higher than for monotonically decreasing flux profiles. 
Note that a locally negative entropy production rate does not contradict the 
second law of thermodynamics as it only states that the entropy of the uni- 
verse should grow, that is, the global — not local — entropy production rate 
should be positive. Nevertheless, imposing a positive local entropy produc- 
tion as an additional constraint in the variational problem, we do recover a 
profile similar to Ozawa and Ohmura (1997). 

Also represented on Fig. [T] are the temperature profile obtained at radia- 
tive equilibrium (that is, all the convective fluxes set to zero) and the con- 
vective flux profile obtained by imposing an isothermal temperature profile. 
These two cases correspond to the bounding profiles for which the entropy 
production rate is zero. Compared to the radiative equilibrium profile, the 
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Figure 2: Potential Temperature Profiles at MEP state (solid lines) and 
at radiative equilibrium (dashed lines) for various longwave optical depths: 
r L = 2 (blue), r L = 3 (red) and tl = 4 (yellow). 

MEP profile has a much more moderate vertical temperature gradient. Sim- 
ilarly, the convective fluxes at MEP state are roughly half those needed to 
maintain a isothermal temperature profile. Note that the temperature pro- 
files so obtain are not statically stable: Fig. [2] shows the vertical profiles of 
(dry) potential temperature, defined by 

, = r(5)", (17) 

for both the radiative equilibrium profiles and the MEP profiles. As ex- 
pected, the radiative equilibrium profiles are not statically stable (potential 
temperature is a decreasing function of height up to roughly 600 hPa). More 
interestingly, the MEP profiles are not statically stable either, as the poten- 
tial temperature is a decreasing function of height up to about 800 hPa - 
which corresponds to the negative entropy production rate zone. Note that 
it could be possible to impose static stability as a constraint in the varia- 
tional problem, possibly with an outcome comparable to imposing a locally 
positive entropy production rate. This constraint would be analogous to the 



method of convective adjustment (Manabe and Strickler 1964). However, we 
shall see in the next section that it is not necessary to do so; replacing the 
gray radiative model by a realistic radiative model turns the MEP profiles 
into statically stable profiles. 

The temperature and convective flux profile obtained at MEP state with 
this simple gray atmosphere radiative model are qualitatively realistic. Nev- 
ertheless, some criticism may still apply. In the first place, the validity of 
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the MEP hypothesis hitherto lies only on empirical evidence. Besides, the 
radiative hypothesis are highly idealized. Diagnostic studies of the entropy 
production rate and other thermodynamic quantities with complex models 
(Pascale et al. 2011) do not fully agree with the MEP hypothesis on the 
vertical dimension. In particular, many different processes contribute to the 
entropy budget (Pauluis and Held 2002b; Pauluis and Held 2002a) and one 
can think of two atmospheres with the same total entropy production rate 
but radically different processes responsible for this entropy production (Volk 
and Pauluis 2010). However, it is not clear how limited the relevance of di- 
agnostic studies with complex models for assessing the validity of the MEP 
hypothesis can be, insofar as the formulation of the variational problem im- 
plies some degree of freedom in the system. In the absence of any unknown 
degree of freedom, that is when all the degrees of freedom are explicitly re- 
solved, possibly through sub-grid parameterizations, there is no room for an 
optimization principle anyway, as the dynamics of the system is fully deter- 
mined by the initial data. Hence, to further inquire on the relevance of the 
MEP variational problem for vertical heat fluxes, we suggest to use an inter- 
mediate approach by solving the variational problem with a realistic radiative 
model, but still no explicit representation of convection. This approach is 
presented in the next section. 



3 MEP Profiles with the NEF model 
3.1 A NEF band model 

In this section, we briefly describe a new radiative model based on the Net 



Exchange Formulation (Dufresne et al. 2005). The purpose of this model is to 



reach a balance between a realistic description of the absorption properties of 
the major radiatively active constituents of the terrestrial atmosphere while 
keeping a relatively smooth dependence of the radiative flux with respect to 
the temperature profile. This last requirement is necessary to use the model 
in the framework of a variational problem. The model presented briefly in 
this section is described in more details in the online supplementary material. 



The fundamental idea of the Net Exchange Formulation (Green 1967) is 
to replace the traditional description of energy exchanges in terms of fluxes 
by a description in terms of rates of energy exchange between two given 
layers. In the flux formulation, the description is local: for a given level in 
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the atmosphere, we are interested in the upward and downward fluxes at 
this precise level. The local laws of emission and absorption allow one to 
write the radiative transfer equation, the solutions of which give the profiles 
of radiative fluxes. On the contrary, in the Net Exchange Formulation, we 
consider two layers in the atmosphere — say i and j — and, using the same 
laws of emission and absorption as usual, we write the net energy exchange 
rate (NER) between these two layers at wavelength u, denoted ip^: 

Va = [ dPi [ dP 3 [ d^a I/ (P i ,j)a v (Pj,j)(B I/ (Tj) — P„(Tj))I„(7), 

(18) 

where r(P i; Pj) denotes the set of all optical paths between point p of layer 
% and point Pj of layer j, Sj the volume of layer % (or a boundary surface - 
the surface of the Earth or space), a u (Pi,j) is the absorption coefficient (1 
if % is a boundary surface rather than a layer), B V (T) the Planck function, 

Iu(l) = exp ^— J k l ,(s)ds\, and k v is the extinction coefficient. The NER 

satisfy simple elementary properties, like antisymmetry: ip- = —ipji, energy 
conservation £\ j ^Xj = (this is a consequence of antisymmetry) or the 
second law of thermodynamics (tp^ has the same sign as Tj — Tj). Here we 

r+oo 

are only interested in the spectrally integrated NER: ipij = / tp^du. 



J% 3 ~ \ 
JO 

In the longwave domain, we decompose the spectrum in 22 narrow bands, 
and in each band we account for absorption by water vapor and carbon 
dioxide only. The absorption coefficient is computed by using the Goody 



(1952) statistical model with the data from Rodgers and Walshaw (1966). 
For the spatial integration, the diffusive approximation is performed with the 
standard diffusion factor /i 1/1.66. Apart from the absorption data, given 
once and for all, the inputs of the model are the water vapor density and 
temperature profiles and carbon dioxide concentration. Of course, one may 
either fix absolute or relative humidity. 

In the shortwave domain, absorption by water vapor and ozone is ac- 



counted for by adapting the parameterization from Lacis and Hansen (1974). 
Apart from the absorption functions determined thanks to laboratory data, 
the input parameters for the model are the water-vapor density and ozone 
density profiles, as well as surface albedo and solar constant. 

The radiative budget of an atmospheric layer is given simply by sum- 
ming over all the NER terms involving the layer in question: for layer i, the 
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radiative budget Ri is given by 

jV+l 

Ri = Y,^i ( 19 ) 

This is quite different from the usual flux formulation for which the radiative 
budget is the divergence of the flux. In particular, it allows for finer discus- 
sions as one may decompose the radiative budget into the NER matrix tp^ 
to see which contributions dominate in the radiative exchanges. 



3.2 The variational problem 

In this section,we consider again a one- dimensional atmospheric column (with 
no lateral fluxes), except that contrary to section [2j the radiative budget is 



now computed using the NEF model presented in section 3.1 The radiative 
budget is no longer expressed as the divergence of a flux as previously, but 
rather directly in terms of the net rate of energy exchanged with all the other 
layers: let us write Ri(T , T\, . . . , T/v) for the radiative budget of layer i and 
7i the convergence of the convective flux (hence 7$ = F c {i — 1) — F c {i)) ■ The 
thermodynamic equation reads: 

dT- 

c p~qI = Ri(T , . . . ,T N ) + v. (20) 

At steady-state, we simply have Ri(T , . . . , T/v) + 7, = in each layer (1 < 
% < N), but 7i is not a known function of Tj. The rate of thermodynamic 
entropy production associated to convection is still given by 

N 

a(T ,...,T N ) = J2^ ( 21 ) 
= _^ R i (T T N ) t (22) 

k=0 1 



with the global constraint 

1 

J]7 fc = $^(To,...,,Tv) = 0. (23) 



N N 



k=0 k=0 



12 



Thus, the MEP variational problem reads 



N 



max 

(To,..., Tat) 



a(T ,...,T N ) I J2 R k(To,...,,T N ) = 0\. (24) 



fc=0 



This variational problem is equivalent to Eq. [16] except that here the radia- 
tive budget is computed with a different radiative model. As a consequence 
it is now more convenient to work with the convergence of convective flux 
rather than the fluxes themselves, which has the consequence that a con- 
straint appears (this constraint is equivalent to Fc{N) = 0). Furthermore, 
writing an equivalent variational problem in terms of the 7j variables would 
imply to invert the radiative budget. Unlike the case of section |2j it is not 
straightforward to do so here, because of a more realistic formulation of the 
radiative exchanges. Note that the advantage of formulating the problem in 
terms of the energy convergence 7 is that it makes clear that the structure 
of the problem is similar to the case of horizontal heat transport (Herbert, 



Paillard, and Dubrulle 2011 Herbert et al. 2011). In particular, the vari- 
ational problem written in this form would remain valid if we considered 
lateral fluxes as well. 



3.3 Results for the NEF model 



Solving the variational problem (24) with standard optimization algorithms, 
we observe that there is a unique solution for a given set of parameters. The 
input parameters for the models are the number of atmospheric layers, the 
water- vapor density profile (fixing either the absolute humidity or the relative 
humidity), the ozone profile, the carbon dioxide concentration and the surface 
albedo (see Table [2j. As the number of vertical layers N increases, the MEP 
temperature and convective flux profiles converge. A good approximation is 
already reached for N = 9 and we shall keep this setting for the rest of this 
paper. For the other parameters, we choose classical values; we use several 



standard water vapor density profiles (McClatchey et al. 1972), and we keep 
the pre-industrial CO2 concentration in this section. 

The temperature and convective flux profiles obtained at MEP state for 
the standard profiles (fixing the absolute humidity) are shown on Fig. |3} 
Compared to the gray atmosphere, we find no inversion in the convection 
profile, that is no region with a locally negative entropy production. The 
vertical temperature gradient is also less steep than for the gray atmosphere. 
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Table 2: Input Parameters for the NEF Radiative Model: 



Parameter 



Symbol 



Value 



Vertical Resolution 
Diffusion factor 

Water vapor profile 
Ozone profile 

Carbon Dioxide concentration 
Solar Constant 
Surface Albedo 



N lev 
H 

PH 2 o{p) 

Po 3 (p) 



PC0 2 

S 

Oin 



9 

1/1.66 
Trop 
SAS 
< SAW 
MLS 
MLW 

v 

280 ppmv 
1368 W.m" 2 
0.1 




240 250 260 270 280 290 20 40 60 80 

Temperature (K) Convective Flux (W.m - -) 



Figure 3: Left: MEP temperature profiles for the standard atmospheric 
profiles: Tropical (blue), Mid-Latitude Summer (red), Mid-Latitude Win- 
ter (yellow), Sub- Arctic Summer (green) and Sub- Arctic Winter (light blue). 
Right: MEP convective flux profiles for the standard atmospheric profiles: 
Tropical (blue), Mid-Latitude Summer (red), Mid-Latitude Winter (yellow), 
Sub- Arctic Summer (green) and Sub- Arctic Winter (light blue). 
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Figure 4: Potential temperature profiles at MEP state for the standard atmo- 
spheric profiles: Tropical (blue), Mid-Latitude Summer (red), Mid-Latitude 
Winter (yellow), Sub- Arctic Summer (green) and Sub- Arctic Winter (light 
blue). 



Note that the vertical temperature gradient is not very sensitive to the water- 
vapor content of the atmosphere: for the various standard profiles, the tem- 
perature profiles are only shifted towards warmer temperatures when the 
water-vapor content increases. Beyond the general structure of the temper- 
ature and convective flux profiles, one may wonder if the MEP profiles are 
statically stable. To this end, we define the (dry) potential temperature as 
in section [2] by 

»-<*)"• < 25 > 

and compute the potential temperature profiles corresponding to the MEP 
profiles discussed above, represented on Fig. |4} In all cases, 6 decreases with 
altitude, which means that the profiles are statically stable. On the contrary, 
the MEP profiles obtained for the gray atmosphere in section [2] were not 
statically stable. It is worth emphasizing that statically stable temperature 
profiles arise here as solutions of the MEP variational problem without spec- 
ifying this as a dynamical constraint. This means that this basic property of 
steady-state profiles is well reproduced by the MEP principle, provided that 
the radiative model is sufficiently realistic. In particular, there is no need 
here to patch the MEP model with a convective adjustment as Pujol and 
Fort ( [2002] ) did for the case of a gray atmosphere. The surface temperature 



obtained at MEP state for the various standard atmospheric profiles is shown 
in Table [3j The results for the tropical and mid-latitude summer profiles are 
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Table 3: Surface temperature and climate sensitivity (surface temperature 
variation for a CO2 doubling experiment) for different water-vapor profiles, 
at MEP state, for the NEF radiative model: 



MEP Surface Temperature MEP Climate Sensitivity 



Tropical 24.8 0.66 

Mid-Latitude Summer 22.4 0.60 

Mid-Latitude Winter 17.1 0.36 

Sub-Arctic Summer 20.8 0.50 

Sub-Arctic Winter 14.6 0.24 

Dufresne and Bony (2008) — 0.4 ± 0.3 



quite reasonable, but the surface temperature for the mid-latitude winter 
and sub-arctic profiles seems relatively high. Two factors can account for 
this high temperature at MEP state: first of all, the effect of clouds is not 
taken into account, and second, the results presented here are all obtained 
with a fixed value for the surface albedo: a g = 0.1. Hence the idealized 
conditions we are using here with the idealized profiles to investigate the 
structure of the MEP profiles is not necessarily meant to represent as closely 
as possible realistic cases. In particular, the arctic profiles are interesting in 
that they feature a lower content of water vapor, rather than because they 
correspond to some geographical location. 

4 Variations in atmospheric composition 

In the previous section, we have shown that replacing the ideal assumptions 
for radiation used in section [2]by a realistic radiative model yields reasonable 
temperature and convective flux profiles. An advantage of this new radiative 
scheme is that it allows for variations of the physical parameters character- 
izing the absorption properties of the atmosphere. Of particular interest are 
the variations of the concentration of carbon dioxide in the atmosphere. As 
a test bed, we will consider in this section the effect on the temperature 
profiles of a doubling of the carbon dioxide concentration. This is a classical 
experiment to which virtually all the climate models are subjected (Randall 
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Figure 5: Temperature response for a CO2 doubling experiment for the stan- 
dard atmospheric (water- vapor) profiles: Tropical (blue), Mid-Latitude Sum- 
mer (red), Mid-Latitude Winter (yellow), Sub- Arctic Summer (green) and 
Sub- Arctic Winter (light blue). 



et al. 2007). One should however remain aware that with the present MEP 
model, very little feedback processes are included, while these feedbacks - 
for instance the water-vapor, albedo, lapse rate or clouds feedbacks — usu- 
ally play a big role in amplifying the climate sensitivity (see e.g. Schlesinger 
and M itchell ( [1987] ; Meehl et al. ( [2004] ); Bony et al. p006| ); Soden and Held 
02006) ). 

Figure [5] shows the response of the temperature profiles obtained with 
the MEP model — using the NEF radiative scheme — for a doubling of the 
CO2 concentration. As expected, in the lower atmosphere, the temperature 
increases when the CO2 concentration increases, while the opposite response 
is observed in the upper atmosphere. The temperature response to CO2 
doubling is greater for profiles with a higher absolute humidity. 

The response of the surface temperature to the variation of the carbon 
dioxide concentration — usually called climate sensitivity — ranges from 
0.24 K (for the sub-arctic winter profile) to 0.66 K (for the tropical profile), 
as shown in table [3] To compare these values with the literature, we need 
to be careful about the feedbacks included in the model we wish to compare 
to. Indeed, if the overall climate sensitivity is still a subject of debate, this is 
mainly due to poorly understood feedbacks, like the cloud feedback (Stephens 



2005), which are not accounted for in the present study. On the contrary, the 
"basic" response of climate to CO2 doubling — due to a change in the Planck 
function (Hansen et al. 1984) — is well agreed-upon: in terms of temperature, 
a rise of 1.2 K is considered an accurate estimate. The "historical" estimate 
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of the seminal study by Manabe and Wetherald (1967) with a radiative model 
with convective adjustment was 1.3 K. Other early studies with radiative- 
convective models with no or little feedbacks provided values as low as 0.8 
K (Hummel and Kuhn 1981b Hummel and Kuhn 1981a) or 0.7 K (Hunt 



1981) — see also Schlesinger and Mitchell (1987). More recently, multi- model 



means with state-of-the-art models separating the contributions of various 
feedbacks (Dufresne and Bony 2008) revealed an estimate of 0.4±0.3 K, when 
taking into account only the effect of the Planck function and the lapse-rate 
feedback, as relevant here. The values obtained with the MEP model for 
idealized atmospheric profiles coincide well with this range. Further studies 
would be necessary to compare more accurately the climate sensitivity of a 
MEP model — possibly including feedbacks like the water-vapor feedback 
(Herbert and Paillard 2012) or the albedo feedback (Herbert, Paillard, and 
Dubrulle 2011) - - to that of radiative-convective models, or even general 
circulation models. 



5 Conclusion 



We have investigated the possibility of computing the vertical transport of 
energy in the atmosphere using a thermodynamic variational principle. The 
variational principle relies on the Maximum Entropy Production conjecture 



(Ozawa et al. 2003; Kleidon and Lorenz 2005; Martyushev and Seleznev 2006 



Dewar et al. 2012), which assumes that the energy fluxes at steady state co- 



incide with the values maximizing the associated thermodynamic entropy 
production. In spite of the lack of theoretical understanding for this hypoth- 
esis, the MEP conjecture has proven useful in several practical applications, 



in particular in climate sciences (Herbert and Paillard 2012). The advantage 



of this formulation is that the energy transport is computed in an implicit 
way, which means that no explicit form is required for the parameterization 
of the sub-grid phenomena effecting the vertical transport in reality. In par- 
ticular, the need for coefficient adjustment is reduced. In addition, it is also 
a very fast solution, numerically speaking. 

In this study, we have tested the realism of the steady-state tempera- 
ture and convective flux profiles obtained by representing the vertical energy 
transport with the MEP principle. Two frameworks have been explored: in 
a first step, we have discussed the highly idealized case of a gray atmosphere 



already considered by Ozawa and Ohmura (1997). In a second step, we have 
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built a model based on realistic assumptions for radiation. Specifically, we 
have designed a radiative scheme based on the Net Exchange Formulation 
(Dufresne et al. |2005[ ) taking into account narrow spectral bands for the ab- 
sorption by water-vapor and carbon dioxide. For several realistic atmospheric 
profiles, we have obtained MEP temperature and convective flux profiles with 
a satisfactory accuracy. We have shown that in this framework, it is possi- 
ble to obtain small but realistic responses for the temperature profiles when 
the composition of the atmosphere is modified. In particular, we have de- 
scribed preliminary results for climate sensitivity obtained in a CO2 doubling 
experiment. 

These results show that in spite of what diagnostic studies based on Gen- 
eral Circulation Models seem to indicate (Pascale et al. 2011), the MEP hy- 
pothesis gives qualitatively realistic profiles when used in conjunction with 
a simple but realistic radiative scheme, even in — and perhaps because of 
- the absence of explicit representation for the hydrological cycle or climate 
feedbacks like clouds, albedo, water- vapor,... Nevertheless, inclusion of such 
processes in a MEP model would be crucial for quantitative applications, and 
in particular for more precise estimations of the climate sensitivity. It also 
paves the way to applications of this type of models to climate problems of 
direct relevance, by varying — at very small computational cost — parameters 
such as carbon dioxide concentration. 



References 

Bony, S., R. Colman, V. Kattsov, R. Allan, C. Bretherton, J.-L. Dufresne, A. 
Hall, S. Hallegatte, M. M. Holland, W. J. Ingram, D. Randall, B. J. Soden, 
G. Tselioudis, and M. J. Webb (2006). "How well do we understand and 
evaluate climate change feedback processes?" J. Climate 19, 3445. 

Dewar, R, C Lineweaver, R Niven, and K Regenauer-Lieb, eds. (2012). Be- 
yond The Second Law: Entropy Production and N on- Equilibrium Systems. 
Understanding Complex Systems. Springer. 

Dewar, R. C. (2003). "Information theory explanation of the fluctuation the- 
orem, maximum entropy production and self-organized criticality in non- 
equilibrium stationary states". J. Phys. A 36, 631-641. 
- (2005). "Maximum entropy production and the fluctuation theorem" . J. 
Phys. A 38, 371-381. 



19 



Dufresne, J.-L. and S. Bony (2008). "An Assessment of the Primary Sources 
of Spread of Global Warming Estimates from Coupled Atmosphere-Ocean 
Models". J. Climate 21.19, 5135-5144. 

Dufresne, J.-L., R. Fournier, C. Hourdin, and F. Hourdin (2005). "Net Ex- 
change Reformulation of Radiative Transfer in the C0 2 15/zm Band on 
Mars." J. Atmos. Sci. 62, 3303-3319. 

Gerard, J.-C, D Delcourt, and L. M. Francois (1990). "The maximum en- 
tropy production principle in climate models: Application to the faint 
young sun paradox" . Q. J. R. Meteorol. Soc. 116, 1123-1132. 

Goody, R. M. (1952). "A statistical model for water vapour absorption". Q. 
J. R. Meteorol. Soc. 78, 165-169. 

Grassl, H. (1981). "The climate at maximum entropy production by merid- 
ional atmospheric and oceanic heat fluxes". Q. J. R. Meteorol. Soc. 107, 
153-166. 

Green, J. S. A. (1967). "Division of radiative streams into internal transfer 
and cooling to space". Q. J. R. Meteorol. Soc. 93, 371-372. 

Grinstein, G and R Linsker (2007). "Comments on a derivation and appli- 
cation of the 'maximum entropy production' principle". J. Phys. A 40, 
9717-9720. 

Groot, S. R. de and P Mazur (1962). N on- equilibrium Thermodynamics. 

North-Holland, Amsterdam. 
Hansen, J., A. Lacis, D Rind, G Russell, P Stone, I Fung, R Ruedy, and 

J Lcrner (1984). Climate Sensitivity: Analysis of Feedback Mechanisms. 

Vol. 29. Geophys. Monogr. Amer. Geophys. Union. 
Herbert, C. and D. Paillard (2012). "Predictive use of the Maximum Entropy 

Production principle for Past and Present Climates". In: Beyond The 

Second Law: Entropy Production and N on- Equilibrium Systems. Ed. by R. 

Dewar, C. Lineweaver, R. Niven, and K. Regenauer-Lieb. Understanding 

Complex Systems. Springer. 
Herbert, C, D. Paillard, and B. Dubrulle (2011). "Entropy production and 

multiple equilibria: the case of the ice-albedo feedback" . Earth Syst. Dy- 

nam. 2, 13-23. 

Herbert, C, D. Paillard, M. Kageyama, and B. Dubrulle (2011). "Present 
and Last Glacial Maximum climates as states of maximum entropy pro- 
duction". Q. J. R. Meteorol. Soc. 137, 1059-1069. 

Hummel, J. R. and W. R. Kuhn (1981a). "An atmospheric radiative-convective 
model with interactive water vapor transport and cloud development". 
Tellus 33, 372-381. 



20 



Hummel, J. R. and W. R. Kuhn (1981b). "Comparison of radiative-convective 
models with constant and pressure-dependent lapse rates". Tellus 33, 
254-261. 

Hunt, B. G. (1981). "An examination of some feedback mechanisms in the 
carbon dioxide climate problem". Tellus 33, 78-88. 

Kleidon, A and R Lorenz, eds. (2005). N on- equilibrium Thermodynamics and 
the Production of Entropy: Life, Earth, and Beyond. Understanding Com- 
plex Systems. Springer, Berlin. 

Lacis, A. A. and J. E. Hansen (1974). "A parameterization for the absorption 
of solar radiation in the Earth's atmosphere". J. Atmos. Sci. 31, 118-133. 

Liou, K. (2002). An Introduction to Atmospheric Radiation. 2nd ed. Interna- 
tional Geophysics. Academic Press. 

Lucarini, V., K. Fraedrich, and F. Ragone (2011). "New Results on the Ther- 
modynamic Properties of the Climate System". J. Atmos. Sci. 68.10, 
2438-2458. 

Manabe, S and R Strickler (1964). "Thermal equilibrium of the atmosphere 
with a convective adjustment". J. Atmos. Sci. 21, 361. 

Manabe, S. and R. Wetherald (1967). "Thermal equilibrium of the atmo- 
sphere with a given distribution of relative humidity" . J. Atmos. Sci. 24, 
241. 

Martyushev, L. M. and V. D. Seleznev (2006). "Maximum entropy production 

principle in physics, chemistry and biology". Phys. Rep. 426, 1-45. 
McClatchey, R., J. Selby, F. Volz, R. Fenn, and J. Garing (1972). Optical 

properties of the atmosphere. Air Force Cambridge Res. Lab. 
Meehl, G. A., W. M. Washington, J. M. Arblaster, and A Hu (2004). "Factors 

affecting climate sensitivity in global coupled models". J. Climate 17.7, 

1584-1596. 

O'Brien, D. M. and G. L. Stephens (1995). "Entropy and climate. II: Simple 

models". Q. J. R. Meteorol. Soc. 121, 1773-1796. 
Onsager, L. (1931). "Reciprocal Relations in Irreversible Processes. I." Phys. 

Rev. 37, 405-426. 

Ozawa, H and A Ohmura (1997). "Thermodynamics of a global-mean state 
of the atmosphere — a state of maximum entropy increase". J. Climate 
10, 441-445. 

Ozawa, H, A Ohmura, R. Lorenz, and T Pujol (2003). "The second law of 
thermodynamics and the global climate system: A review of the maximum 
entropy production principle". Rev. Geophys. 41, 1018. 



21 



Paltridge, G. W. (1975). "Global dynamics and climate-a system of minimum 
entropy exchange" . Q. J. R. Meteorol. Soc. 101, 475-484. 

Pascale, S., J. M. Gregory, M. H. P. Ambaum, and R. Tailleux (2011). "A 
parametric sensitivity study of entropy production and kinetic energy 
dissipation using the FAMOUS AOGCM". Clim. Dyn. 38, 1211-1227. 

Pascale, S., J. M. Gregory, M. H. P. Ambaum, R. Tailleux, and V. Lucarini 
(2012). "MEP solution for a minimal climate model: success and limita- 
tion of a variational problem". Earth Syst. Dynam. 3, 19-32. 

Pauluis, O. (2011). "Water vapor and mechanical work: a comparison of 
Carnot and steam cycles." J. Atmos. Sci. 68, 91-102. 

Pauluis, O. and I. M. Held (2002a). "Entropy budget of an atmosphere in 
radiative-convective equilibrium. Part I: Maximum work and frictional 
dissipation". J. Atmos. Sci. 59.2, 125-139. 
- (2002b). "Entropy budget of an atmosphere in radiative-convective equi- 
librium. Part II: Latent heat transport and moist processes". J. Atmos. 
Sci. 59.2, 140-149. 

Pujol, T. (2003). "Eddy heat diffusivity at maximum dissipation in a radiative- 
convective one-dimensional climate model". J. Meteor. Soc. Japan 81.2, 
305-315. 

Pujol, T. and J. Fort (2002). "States of maximum entropy production in a 
one-dimensional vertical model with convective adjustment" . Tellus A 54, 
363-369. 

Randall, D. A., R. A. Wood, S Bony, R Colman, T Fichefet, J Fyfe, V 
Kattsov, A Pitman, J Shukla, J Srinivasan, R. J. Stouffer, A Sumi, and 
K. E. Taylor (2007). "Climate Models and Their Evaluation". In: Cli- 
mate Change 2007: The Physical Science Basis. Contribution of Working 
Group I to the Fourth Assessment Report of the Intergovernmental Panel 
on Climate Change. Ed. by S Solomon, D Qin, M Manning, Z Chen, M 
Marquis, K. B. Averyt, M Tignor, and H. L. Miller. Cambridge University 
Press, Cambridge, United Kingdom and New York, NY, USA. 

Rodgers, C. D. and C. D. Walshaw (1966). "The computation of infra-red 
cooling rate in planetary atmospheres" . Q. J. R. Meteorol. Soc. 92, 57-92. 

Schlesinger, M. E. and J. F. B. Mitchell (1987). "Climate model simulations 
of the equilibrium climatic response to increased carbon dioxide". Rev. 
Geophys. 25, 760-798. 

Soden, B. J. and I. M. Held (2006). "An assessment of climate feedbacks in 
coupled ocean-atmosphere models". J. Climate 19, 3354-3360. 



22 



Stephens, G. L. (2005). "Cloud feedbacks in the climate system: A critical 
review" . J. Climate 18, 237-273. 

Stone, P. H. (1978). "Constraints on dynamical transports of energy on a 
spherical planet". Dyn. Atmos. Oceans 2, 123-139. 

Volk, T and O. Pauluis (2010). "It is not the entropy you produce, rather, 
how you produce it". Phil. Trans. R. Soc. B 365, 1317-1322. 

Wang, B., T. Nakajima, and G. Shi (2008). "Cloud and Water Vapor Feed- 
backs in a Vertical Energy-Balance Model with Maximum Entropy Pro- 
duction". J. Climate 21.24, 6689-6697. 



23 



